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ABSTRACT 

We have investigated the pulsar PSR B2224+65 and its X-ray jet with XMM-Newton. Apart 
from the long X-ray jet which is almost perpendicular to the direction of proper motion, a 
putative extended feature at the pulsar position, which oriented in the opposite direction of the 
proper motion, is also suggested by this deep X-ray imaging. Non-detection of any coherent 
X-ray pulsation disfavors the magnetospheric origin of the X-rays observed from the position of 
PSR B2224-I-65 and hence suggest that the interpretation of pulsar wind nebula is more viable. 
We have also probed the origin of PSR B2224-I-65 and identified a runaway star, which possibly 
originated from the Cygnus 0B9 assoication, as a candidate for the former binary companion of 
the neutron star's progenitor. 
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INTRODUCTION 
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While the signature of rotation-powered pul- 
sars is their pulsed emission, bulk of the rota- 
tional energy is going to power a relativistic out- 
flow of magnetohydrodynamic wind particles with 
Lorentz factors up to 7 ^ 10^. Although pulsar 
winds cannot be observed directly, their presence 
can be inferred through their interaction with the 
interstellar surroundings. If the wind reaches a 
termination shock, a pulsar wind nebula (PWN) 
is formed. The term PWN had been originally 
referred to the shocked emission associated with 
young and powerful pulsars (e.g. Crab) where the 
wind is effectively confined by the magnetic fields 
of the surrounding supernova remnants and rela- 
tivistic leptons reach high energy densities which 
make these nebulae usually luminous synchrotron 
sources in the radio and X-ray regimes. 

A different class of PWNe which is associated 
with older and less powerful pulsars is bow-shock 
nebulae. These nebulae are not embedded in su- 
pernova remnants but interact directly with the 
interstellar medium (ISM) . If the pulsar has a high 
velocity, the nebula is confined by the ram pres- 
sure of the ISM and forms the typical bow-shock 
structure. The bow-shock size is determined by 



the balance between the wind and the ISM ram 
pressure. Even in this old and less powerful pul- 
sars, there is a significant fraction of the whole 
pulsar spin-down power which is going to feed the 
pulsar wind (see Bucciantini & Bandiera 2001). 

Among all the pulsars with PWNe detected, 
PSR B2224+65 is one of the most interesting sys- 
tem. This pulsar has a rotational period and 
a period derivative of P ~ 0.68 s and P 
9.6 X 10""'^^ s s'-^, respectively. This suggests 
PSR B2224-h65 has a spin-down age of 1.13 x 10^ 
years. The spin parameters imply a spin-down 
luminosity and a surface dipole magnetic field 
strength oi E ^ 1.2 x 10'^'^ erg s^^ and B ~ 
2.6 x 10^^ G respectively. Concerning its distance, 
while the radio dispersion measure implies a value 
of 1.86 kpc (cf. Cordes & Lazio 2002), Chatterjee 
& Cordes (2004) found that by modeling the head 
of the shocked emission in Ha the pulsar distance 
could be as close as 1 kpc. 

Hui & Becker (2007) have examined its X-ray 
properties with the Chandra data taken in Octo- 
ber 2000. Brief results of the same observation 
have also been reported by Wong et al. (2003) 
and Zavlin & Pavlov (2004). Despite the pho- 
ton statistics is limited in this data set, the spec- 
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tral analysis suggests that the bulk of the X-ray 
emission from PSR B2224-I-65 is non-thermal. An 
independent analysis that includes also the sec- 
ond epoch Chandra observations taken in August- 
October 2006 has confirmed the non-thermal na- 
ture of the X-rays detected at the pulsar position 
(Johnson & Wang 2010). Nevertheless, the large 
frame time (i.e. 3.2 sec) of these Chandra obser- 
vations do not allow us to constrain the temporal 
properties of the pulsar. Therefore, with the ex- 
isting data, we cannot determine whether the bulk 
of these observed X-rays are originated from the 
pulsar magnetosphere or from the PWN. 

The proper motion of PSR B2224-I-65 is 
among the largest in the currently known pul- 
sar population with a transverse velocity of 
V = 865(cZ/l kpc) km s~^ (Harrison, Lyne & An- 
derson 1993). Its shock interaction with the ISM 
produces the most well-known Guitar Nebula as 
observed in Hq (Cordes et al. 1993). Although 
no evidence for the X-ray emission can be found 
in the region of the Ha nebula, an X-ray jet ap- 
parently associated with PSR B2224-I-65 has been 
detected by Chandra (Hui & Becker 2007; Wong 
et al. 2003). The most peculiar fact of this jet 
is that its orientation deviates by ~ 118° from 
the pulsar's proper motion direction. Its X-ray 
spectrum is also found to be non-thermal. 

Kargaltsev & Pavlov (2008) have argued that 
this jet-like feature may be associated with a 
nearby point source instead of with PSR B2224+65 
However, this possibility has been ruled out by a 
recent analysis of the multi-epoch Chandra obser- 
vations. Johnson & Wang (2010) have found that 
the sharp leading edge have a proper motion con- 
sistent with that of the pulsar, which secures the 
association between the jet and PSR B2224-I-65. 
Different pictures have been proposed to explain 
this X-ray jet (e.g. Bandiera 2008; Johnson & 
Wang 2010). However, there is still no consensus 
on its nature. 

To further investigate the nature of the X-rays 
from the position of PSR B2224-I-65 and its enig- 
matic jet, a dedicated follow-up observation with 
improved temporal resolution and photon statis- 
tic is needed. In this paper, we present a deep 
XMM-Newton observation in order to study this 
extreme system in details. In §2, we describe the 
details of the XMM-Newton observation and the 
data reduction. In §3 and §4, we present the re- 



sults of the X-ray analysis and the search for the 
7— ray emission from this system respectively. Fi- 
nally, we discuss the physical implications of the 
observed results and search for the birth place of 
this fast-moving neutron star in §5. 

2. X-RAY OBSERVATION WITH XMM- 
NEWTON 

PSR B2224-F65 was observed by XMM-Newton 
on 13-14 July 2009 (Observation ID: 0604420101) 
with an exposure time of 75 ks. The mean epoch 
of this observation is at MJD 55025. The data 
used in this investigation were obtained with 
the European Photon Imaging Camera (EPIC) 
aboard XMM-Newton (Jansen et al. 2001). EPIC 
consists of two Metal Oxide Semiconductor (MOSl/2) 
CCD detectors (Turner et al. 2001) of which half 
of the beam from two of the three X-ray tele- 
scopes is reflected to. The other two halves of the 
incoming photon beams are reflected to a grating 
spectrometer (RGS) (den Herder et al. 2001). The 
third of the three X-ray telescopes is dedicated to 
expose the PN CCD detector solely (Striider et 
al. 2001). The PN CCD was operated in small- 
window mode with a thin fllter to block optical 
stray light. All the events recorded by PN cam- 
era are time tagged with a temporal resolution of 
5.7 ms, which enable us to constrain the temporal 
properties of this pulsar in X-ray regime for the 
flrst time. On the other hand, the MOSl/2 CCDs 
were setup to operate in full-window mode with a 
thin filter in each camera. 

The aimpoint of the satellite in this observa- 
tion is RA=22'>25™52.43'' and Dec=-F65°35'34.0" 
(J2000). With the most updated instrumental 
calibration, we generate the event lists from the 
raw data obtained from all EPIC instruments with 
the tasks emproc and epproc version 9.0.0 of the 
XMM Science Analysis Software (XMMSAS ver- 
sion 9.0.0). Examining the raw data from the 
PN CCD, we did not find any timing anomaly 
which was observed in many of the XMM-Newton 
data sets (cf. Hui & Becker 2006 and references 
therein). This provides us with opportunities for 
an accurate timing analysis. We then created fil- 
tered event files for the energy range 0.5 keV to 
10 keV for all EPIC instruments and selected only 
those events for which the pattern was between 
0-12 for MOS cameras and - 4 for the PN 
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camera. We further cleaned the data by accepting 
only the good times when sky background was low 
and removed all events potentially contaminated 
by bad pixels. After the filtering, the effective ex- 
posures are found to be 46 ks for each MOS camera 
and 32 ks for PN. 

3. ANALYSIS OF XMM-NEWTON OB- 
SERVATION 

3.1. Spatial Analysis 

X-ray images of the 6' x 6' field around 
PSR B2224-K65 as obtained by MOSl/2 and PN 
camera are displayed in Figure [TJ Both image 
have been firstly produced with a binsize of 1" 
which is subsequently smoothed with a Gaussian 
with a kernel size of cr = 3". Both the pulsar and 
the jet can be clearly seen. 

For investigating the spatial variation of the X- 
ray intensity of the jet, we computed the bright- 
ness profiles along its orientation with MOSl/2 
and PN data which are shown in Figure [2] We es- 
timate the counts in the consecutive boxes with 
a size of 25" x 10" from the raw images with 
a binsize of l" along the jet. The orientation 
of these sampling regions are shown in the in- 
sets of Figure [H To estimate the background 
level, we have sampled the source-free regions 
around PSR B2224+65 within a 6' x 6' field-of- 
view. The average background level and its la 
uncertainties are indicated by the horizontal solid 
line and dotted lines in Figure [2] respectively. Ap- 
parently, the jet extends for ~ 2.5' before it falls 
to the estimated background level. Apart from 
the peak corresponds to PSR B2224-f 65, another 
peak at ~ 70" from the pulsar position is also 
observed in both MOSl/2 and PN brightness pro- 
files. Hereafter, we refer this as a "knot-like" fea- 
ture. To determine the position of the knot-like 
feature on the jet, we have performed source de- 
tection individually on each EPIC data set with 
the aid of the XMMSAS task edetect_chain. The 
mean X-ray position inferred from this observation 
is RA=22i>25'"42.24" Dec=-h65°36'00.0" (J2000) 
with a resultant uncertainty of 5 by combining 
the statistical errors inferred from each camera in 
quadrature. 

To further constrain the nature of the knot- 
like feature, we make use of the archival data 
obtained by the Hubble Space Telescope (HST). 



Utilizing the DAOphot catalog (version 4.0) re- 
leased on the Hubble Legacy Archiv^ we have 
identified a possible optical counterpart with a 
magnitude of i? - 22.31 at RA=22'^25'^42.14" 
Dec=-h65°36'01.2" (J2000) which locates within 
the X-ray positional uncertainty of the knolH. The 
root mean square (rms) positional uncertainty is 
^ 0.2 of the optical source in each coordinate 
based on the astrometric correction with 53 refer- 
ence sources in Guide Star Gatalog 2 (Lasker et 
al. 2008). Assuming the optical source is associ- 
ated with the X-ray knot, we can probe its possible 
nature together with X-ray spectral analysis (see 
§3.2). 

With the archival data obtained by Ghan- 
dra X-ray Observatory, it is possible for us to 
place a constraint on the knot's proper motion. 
We utilized the Chandra Interactive Analysis 
of Observations software (GIAO 4.3) for the re- 
analysis. We have reprocessed all the available 
dataof PSRB2224+65 (Obs. IDs 755; 6691; 7400) 
obtained by the Advanced GGD Imaging Spec- 
trometer (AGIS) onboard Ghandra with the most 
updated calibration GALDB (ver. 4.4.5). Differ- 
ent from the analysis reported by Johnson & Wang 
(2010) and Hui & Becker (2007), we have also 
applied the subpixel event repositioning during 
reprocessing the data in order to improve the po- 
sitional accuracy of each event (cf. Li et al. 2004). 
For this analysis, we considered only the Obs. IDs 
755 and 7400 which have an epoch separation of 
~ 6 yrs. The knot can be detected in both images 
with a wavelet source detection algorithm. Since 
Johnson & Wang (2010) found that the rotation 
between these images is negligible, simple transla- 
tional astrometric corrections have been applied. 
The proper motion of the knot inferred from these 
two frames is /ira = —1.0 ± 0.4 arcsec/yrs and 
MRA = 0.18 ± 0.15 arcsec/yrs. These errors (Icr) 
are computed by combining the X-ray positional 
errors of the knot in both frames and the corre- 
sponding astrometric uncertainties in quadrature. 
The error budget is dominated by the uncertain- 
ties in determining the positions of the knot. In 
view of the large errors, the proper motion of the 
knot remains unconstrained. 



^http:/ /hla. stsci.edu/hlaview. html 

^The position of this optical counterpart is inferred by a 
HST observation performed on 2001-11-22 (Proposal ID; 
9129) 
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The insets of Figure [T] show the 1' x 1' close- 
up view of PSR B2224-I-65 in the corresponding 
cameras. The X-ray emission from the pulsar po- 
sition appears to be slightly elongated for ^ 8" 
behind proper motion direction as illustrated by 
the arrow in Figure [T] In both MOSl/2 and FN 
images, this putative feature apparently has the 
same orientation. However, the elongation of this 
feature is close to the edge of the spatial resolution 
of XMM-Newton. This has motivated us to reex- 
amine the archival Chandra data which can offer 
an improved angular resolution. We produced a 
merged epoch X-ray image after correcting the rel- 
ative astrometric offsets among each data set. The 
epoch-merged raw X-ray image of FSR B2224-I-65 
is shown in Figure [3^. Figure [3}d shows the X-ray 
brightness profile as measured by eight consecu- 
tive boxes with the orientation illustrated in the 
inset. Interestingly, both image and the bright- 
ness profile suggest a slight deviation from a point 
source with deformation toward a similar direc- 
tion as observed by XMM-Newton. The extent of 
^ 3" is observed in the ACIS data, which can pos- 
sibly be the part of the highest surface brightness, 
though it is not conclusive with the current data 
(see §5 for a further discussion). 

3.2. Spectral Analysis 

With the aid of the XMMSAS task epatplot, 
all the EPIC data are found to be not affected 
by CCD pileup. For PSR B2224+65, we ex- 
tracted its source spectrum from circles with a 
radius of 20" centered at the point source posi- 
tion in both MOSl/2 and FN cameras respec- 
tively. This choice of extraction regions corre- 
sponds to the encircled energy fraction of ~ 75% 
and ~ 80% in MOSl/2 and FN respectively S 
For the jet, we sampled its source spectrum from 
a 122" X 36" box centered at RA=22h25'"40.69" 
Dec=-h65°36' 07.64" (J2000) and the longer side 
oriented 66° east from the north so as to cover 
the whole jet. The background spectra were ex- 
tracted from the nearby regions which are source- 
free and with sufficient size to enable a less bi- 
ased sampling. Response files were computed for 
all datasets by using the XMMSAS tasks rmfgen 
and arfgen. After the background subtraction, we 



http://xmm.esac.esa.int/external/xmm_ user_ sup- 
port/ documentation / uhb /nodol 7. html 



have 74 ± 11 cts (MOSl), 67 ± 10 cts (M0S2), 
189 ± 17 cts (FN) collected from FSR B2224+65 
and 183 ± 18 cts (MOSl), 212 ± 18 cts (M0S2), 
453 ± 30 cts (FN) collected from the jet respec- 
tively. 

We have firstly examined the data obtained 
from each camera individually and found that 
all the spectral parameters inferred from differ- 
ent cameras are consistent within la uncertainties. 
In order to tightly constrain the spectral parame- 
ters, we fitted the data obtained from three cam- 
eras simultaneously for the individual analysis of 
FSR B2224+65 and its jet. According to the pho- 
ton statistics, we grouped each spectrum dynami- 
cally so as to achieve the same signal-to-noise rate 
in each dataset. The adopted spectral binning en- 
sures a sufficiently small errors in each energy bin 
for better discriminating the competing scenarios. 
All the spectral fits were performed with XSFEC 
12.6.0 in the energy range of 0.5 — 10 keV. All the 
quoted errors are la for 1 parameter of interest 
(i.e. Ax^ = 1.0 above the minimum). 

For FSR B2224+65, we found that a simple 
absorbed power-law model yields a goodness-of- 
fit of — 3.35 with 8 d.o.f.. Examining the 
power-law with different spectral binning factor, 
the model remains to provide a very good descrip- 
tion of the data. The best- fit power- law model and 
the residuals for the pulsar spectrum are shown 
in Figure 2] This model yields a column den- 
sity of uh = (1.31^2) ^ 10^^ cm~2, a photon 
index of Fx = 2.0^g'4 and a normalization at 1 
keV of (6.7j:^:g) X 10"^ photons keV'^ cm'^ s"!. 
We have also examined whether the X-rays from 
FSR B2224-I-65 can be described in a thermal- 
dominant emission scenario by fitting the spec- 
trum with a blackbody model. The best-fit model 
yields a column density of nn < 4.6 x 10^° cm~^, 
a temperature of kT = 0.5 ±0.1 keV and an emit- 
ting radius of i?bb = (14-lt2;3) ^ '^ikpc m, where 
c^ikpc at the unit of 1 kpc. The inferred emit- 
ting region is smaller than the radius of the polar 
cap (i.e. Rpc = 175 m) as defined by the last 
open field lines of a dipolar magnetic field. Also, 
the goodness-of-fit resulted from the best-fit black- 
body model is somewhat worse (x^ = 9.28 with 8 
d.o.f.) in comparison with the power-law fit. Ex- 
amining the fitting residuals, some systematic de- 
viations are found between the observed spectrum 
and the best-fit blackbody model. In view of these. 
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the blackbody model is not favored. Therefore, 
we conclude that the simple absorbed power-law 
model is favored in modeling the observed spec- 
trum of PSR B2224-I-65, which agrees with the 
inference drawn by Hui & Becker (2007) and John- 
son & Wang (2010). 

For the jet, we found that a power-law model 
yields a column density of uh = (2.7lQg) x 
10^^ cm~^, a photon index of Tx = 1.2 ± 0.1 and 
a normalization at 1 keV of (Q.ll^'^) x 10^^ pho- 
tons keV"""^ cm~^ s~^. The model results in a 
goodness-of-fit of = 38.37 with 40 d.o.f. The 
inferred column density from fitting the jet spec- 
trum is consistent with that inferred from the pul- 
sar spectrum. On the other hand, in examining 
whether a thermal scenario is capable to describe 
the observed data, we attempted to fit a thermal 
bremsstrahlung model to the spectrum of the jet. 
However, we found that it results in an unphysi- 
cally high temperature of kT ^ 180 keV. There- 
fore, our investigation further supports a non- 
thermal origin of the X-rays from the jet. 

In order to constrain the spectral properties of 
this system more tightly, we follow the method 
adopted by Johnson & Wang (2010) to jointly 
fit individual power-law models for the X-ray 
spectra of the pulsar and the jet with the col- 
umn density in the individual model tied to- 
gether. The results are summarized in Table [T] 
The joint analysis yields a column density of 
nn = (2.51J;?) x lO^i cm'^. The photon index 
and the normalization at 1 keV for the pulsar are 
found to be Tx 2.2l° j and (8.3l ^ g) x 10^^ pho- 
tons keV~^ cm^^ s~^ respectively. The corre- 
sponding best-fit spectral parameters for the jet 
are Tx = 1.2 ± 0.1 and (8.7^};^) x 10"^ pho- 
tons keV^^ cm^^ s~^ respectively. In the en- 
ergy band of 0.5 — 10 keV, the observed flux 
and absorption-correct fluxes for PSR B2224+65 
are found to be (2.4l:^;o) x lO^^-* ergs cm^^ g-i 
and (3.4ti[;Q) x 10~^* ergs cm"^ s~^ respec- 
tively. And the observed flux and absorption- 
corrected flux for its jet in the same energy 
band are (9.112^) 10^^^ ergs cm~^ s~^ and 
(10.0^2 7) ^ 10"^'' ergs cm~^ s~-^ respectively. 

X-ray emission of the jet can possibly be pow- 
ered by the energy loss of the synchrotron emitting 
leptons during their rides from the pulsar. If the 
synchrotron cooling time is less than the timescale 
for the flow of leptons across the feature, steepen- 



ing of the photon index for the synchrotron radia- 
tion is expected along the jet. In order to investi- 
gate this possible spectral variation along the jet, 
we divide the extraction region for the whole jet 
equally into two parts along its orientation. In the 
following, we refer the segment closer to the pulsar 
as region 1 and the further one as region 2. The to- 
tal photon statistics collected by three EPIC cam- 
eras for region 1 and region 2 are 379 ± 28 cts and 
376 ± 26 cts respectively. For modeling the spectra 
of these segments, we also consider the power-law 
model. For a more constraining analysis, we again 
jointly fit the spectra of regions 1 and 2 as well 
as the whole jet with the column density in the 
individual model tied together. The results are 
summarized in Table [T] In this joint analysis, the 
X-ray spectra in both regions can be well-modeled 
with a power-law with a slope Tx '^1.3 and hence 
no variation of photon index have been found. 

Since a knot-like X-ray feature along the jet is 
identified and a possible optical source is found 
within its positional uncertainty, this leads us to 
speculate whether this is a field object just hap- 
pens to locate in the error circle by chance instead 
of intrinsically related to the jet. Therefore, we 
also attempt to examine the X-ray spectrum of 
this knot-like feature. Within a circular region 
with a radius of 15" centered at its nominal posi- 
tion (see §3.1), there are 100 ± 15 cts collected by 
all three EPIC cameras after background subtrac- 
tion. We have examined its spectrum by fitting 
with various single component models which in- 
clude a power-law, a thermal bremsstrahlung as 
well as a collisional ionization equilibrium plasma 
model. In view of the limited photon statistic of 
this feature, we simply assume the column absorp- 
tion is at the level constrained in the aforemen- 
tioned joint analysis (i.e. uh = 2.8 x 10^^ cm~^). 
We found that all three tested models result in a 
similar goodness-of-fit {x^ = 18.1 with 23 d.o.f.) 
and hence we are not able to discriminate these 
competing scenarios. We notice that all the mod- 
els yield an observed flux at the level of fx ^ 
1.5 X 10-14 ergs cm'^ s^^ (0.5 - 10 keV). In 
the HST image, the possible optical counterpart 
of this feature has a magnitude of i? ~ 22.31. 
This gives rise to an X-ray-to-optical flux ratio of 
fx/ Jr 0.77. This value is higher than that ex- 
pected for a field star which typically has a ratio 
of fx/fopt < 0.3 (Maccacaro et al. 1988) but con- 
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forms with the expected range of an AGN which 
typically has a ratio of fx/fopt < 50 (Stocke et al. 
1991). 

Since the nature of the knot-like feature cannot 
be confirmed, we have also reexamined the spec- 
tral behavior of the jet with the knot removed as 
if it were a background source. In particular, as 
the knot is located in region 1, we would like to in- 
vestigate whether its removal has an effect in the 
spectral shape of this region. After subtracting 
the contribution from the knot, we re-do the joint 
spectral fitting of the spectra of regions 1 and 2 as 
well as the whole jet. The results are also summa- 
rized in Table [TJ The photon indices in these two 
regions remain to be consistent within la uncer- 
tainties. Therefore, we conclude that no evidence 
for spectral steepening along the X-ray jet can be 
found in this observation. 

3.3. Temporal Analysis 

3.3.1. Search for the coherent pulsation from 
PSR B 2224 +65 

In order to search for X-ray pulsations from the 
pulsar PSR B2224-I-65, only the data taken from 
the PN small window mode with a temporal res- 
olution of 5.7 ms can be used. We extracted 1450 
counts from a circle with a radius of 20" centered 
at the radio timing position. We further selected 
these events within the good time intervals with 
the counts rate of the entire field lower than 0.02 
cts/sec in order to reduce the contamination from 
the sky background. After filtering, a total of 523 
source counts was left in our study. The pho- 
ton arrival times were corrected to the solar sys- 
tem barycenter with the BARYCEN tool (version: 
1.18; JPL DE200 Earth ephemeris) of XMMSAS. 

We noticed that the pulsar has a glitch at 
MJD ~54266, which is about two years before 
our observation (MJD 55025). Since there is no 
further information about the relaxation time, we 
then performed a detailed search for a range that 
cover both the post-glitch frequency as well as 
the extrapolated value at the mean epoch of the 
XMM-Newton observation. We applied the 
statistics (Buccheri et al. 1983) with the harmon- 
ics number (n) from one to ten and the H-Test (de 
Jager, Raubenheimer & Swanepoel 1989) with the 
searching step of 1/5 corresponding Fourier width 
in the range of 1.4650 to 1.4652 s'^ No significant 



signal was detected. We also folded the light curve 
by using the radio spin period extrapolated for the 
epoch of the XMM-Newton observation. However, 
no meaningful light curve was obtained. The pul- 
sar ephemeris reported from the ATNF Catalogue 
(Manchester et al. 2005), / = 1.46511023680 Hz 
and / = -2.0737 x lO^^^^ sec-^ (at MJD = 
54420.0) were used in our study. We further com- 
puted the upper limit of the pulsed fraction fp 
by fp = p{Nt-Nb)/Nt, where Nt and Nb are 
the number of total photons and the background 
photons respectively and /S is the duty cycle. As- 
suming /3 = 0.5, we placed a la upper Hmit for 
the pulsed fraction to be ~ 8%. 

3.3.2. Search for long-term variability from the 
jet 

As an X-ray jet is possibly subjected to mag- 
netohydrodynamic instability, variability has been 
seen in many PWN systems (e.g. DeLaney et al. 
2006). In order to investigate the long-term vari- 
ability from the jet, we analyse the X-ray obser- 
vations of this field in different epochs. We have 
ignored one of the Chandra observations (Obs ID. 
6691) as its short exposure (~ 10 ks) does not 
provide sufficient photon statistic for this anal- 
ysis. The multi-epoch spectral results are sum- 
marized in Table [21 Since the flux estimates are 
sensitively to the column absorption, we fixed 
the riH in all observations at the value inferred 
from the analysis of the XMM-Newton observa- 
tion as it provides a superior photon statistic in 
the soft band and hence put a better estimate on 
the hh- Also, for the sake comparison, we com- 
pute the flux in an energy range of 0.5 — 10 keV 
for all cases. We found that both observed and 
absorption-corrected fluxes, as well as the photon 
index are consistent among the observations in dif- 
ferent epochs within la uncertainties. Therefore, 
we do not find any evidence of the long-term spec- 
tral and flux variability. 

Apart from searching for the long-term variabil- 
ity, we have also investigated if there is any no- 
ticeable flux variation within each individual ex- 
posure. By examining the light curves, we also 
do not flnd any evidence for the variability within 
each observation window. 
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4. SEARCH FOR 7 RAY EMISSION 
WITH FERMI LARGE AREA TELE- 
SCOPE 

PSR B2224+65 has a spin-down age of 
10^ yrs. It has been speculated that the accel- 
erating regions can still be sustained in the mag- 
netospheres of some old pulsars (e.g. Zhang et 
al. 2004; Hui & Becker 2007). Indeed, pulsed 
7— rays with energies > 0.1 GeV have recently 
been detected from an old radio-quiet pulsar, 
PSR J1836-I-5925, that has a characteristic age 
of 1.8 million years (Abdo et al. 2010). Motivated 
by this discovery, we also attempt to search for 
the 7— ray emission from PSR B2224+65. 

The continuous release of the photon data ob- 
tain in the all-sky survey with the Large Area 
Telescope (LAT) on board the Fermi Gamma-ray 
Space Telescope enables us to search for the 7— ray 
emission from PSR B2224-K65 in MeV-GeV 
regime. In this analysis, we used the data ob- 
tained between 2008 August 4 and 2011 Oct 2. 
We used the Fermi Science Tools v9r23pl package 
to reduce and analyze the data in the vicinity of 
PSR B2224-I-65. Only events that are classified as 
the "Source" clas^ were used. To reduce the con- 
tamination from Earth albedo 7-rays, we excluded 
events with zenith angles greater than 100°. The 
instrument response functions "P7SOURCE_V6" 
were used. Events with energies between 100 MeV 
and 20 GeV were used for our analysis. 

With the aid of the task gtlike, we carried 
out an unbinned maximum-likelihood analysis for 
a region-of-interest (ROI) with a 5° radius cen- 
tered on the position of PSR B2224-I-65 in the en- 
ergy range of 0.1 — 20 GeV. We subtracted the 
background contribution by including the Galac- 
tic diffuse model (gaL2yearp7v6_v0.fit) and the 
isotropic background (iso_p7v6source.txt), as well 
as all sources in the second Fermi/LAT cata- 
log (2FGL; Abdo et al. 2011) within 10° from 
PSR B2224-I-65. We assumed a power-law spec- 
trum for PSR B2224+65 as well as aU 2FGL 
sources in our consideration. Making use of the 
model parameters resulted from the maximum- 
likelihood analysis, we computed the test-statistic 
(TS) map of 2° x 2° centered at the nominal 



An event class in Pass 7 (i.e. evclass=2) that provides 
the most suitable photon data for point source spectral 
analysis. 



position of PSR B2224-h65. The TS valucEl at 
the position of PSR B2224-h65 is found to be 
~ 40. However, we speculated that this can 
possibly be contributed by two previously un- 
known nearby 7— ray sources with their TS val- 
ues of 29 and ~ 38 peaked at RA=336.82° 
Dec=65.95° (J2000) and RA=335.25° Dec=64.98° 
(J2000) respectively. Appending these two ad- 
ditional sources in the model, we have repeated 
the maximum-likelihood analysis and the TS map 
computation. With the contribution from these 
two serendipitous sources subtracted, the TS value 
at the position of PSR B2224-I-65 essentially drops 
to zero and no residual emission can be identified 
in the TS map. Therefore, we conclude that we 
do not find any evidence for 7— ray emission from 
PSR B2224-I-65 in this study. Adopting a power- 
law photon index of F = 2, we place a Icr upper- 
limit of - 4 X 10^12 gi-g cni-2 in 0.1 - 20 GeV. 

5. DISCUSSION 

5.1. Nature of the X-ray emission from 
PSR B2224+65 & its jet 

We have presented a detailed investigation of 
the X-ray emission from PSR B2224-f65 and its 
jet with XMM-Newton. The non-detection of any 
coherent X-ray pulsation suggests that the mag- 
netospheric origin is questionable. In the con- 
text of outer-gap model, the non-thermal pulsed 
X-rays are resulted from the back-flowing charge 
particles from the outer-gap (e.g. Cheng & Zhang 
1999). Based on this model, the electromagnetic 
cascade developed in this flow should give rise to 
an X-ray photon index < 2. The photon index 
of the X-rays from PSR B2224-f65 is 2.2^°;^. The 
best-fit value is found to larger than that expected 
from the magnetosplieric model, though they can 
be reconciled within la uncertainties. 

On the other hand, with the deep X-ray imag- 
ing of PSR B2224-I-65, we have revealed an in- 
teresting feature at the pulsar position which is 
elongated along direction of proper motion as ob- 
served in different cameras (cf. insets in Fig. [ij. 
Although we cannot confirm its extent unambigu- 
ously with the existing data, its orientation is in- 
teresting and worth discussing. X-ray tails have 



^TS value corresponds to the detection significance of ct Ri 
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been found to be associated with many pulsars 
(see Kargaltsev & Pavlov 2010 for a recent review). 
The tails are usually oriented in the opposite direc- 
tion of pulsar motion and behind the bow-shock, 
therefore it is interpreted as the flow of parti- 
cles coming out of the ram-pressure dominated 
PWN. In this scenario, the X-ray emission from 
PSR B2224-I-65 is originated from the synchrotron 
radiation from the wind particles accelerated at 
the termination shock. This is consistent with the 
unpulsed non-thermal emission inferred from this 
observation. For the synchrotron emitting leptons 
of Lorentz factor 7 distribute as N (7) oc 7"'', the 
observed photon index of F ~ 2.2 from the pul- 
sar location is consistent with the value expected 
from the fast-cooling scenario, F = (p + 2) /2, for 
the electron index lies within a typical range of 
p ~ 2 - 3 (Cheng, Taam, & Wang 2006). 

This putative feature, which has an extent 
of ~ 8" in EPIC images, should be resolved 
by Chandra. Nevertheless, in reexamining the 
archival Chandra data, its extent cannot confirm. 
On the other hand, the photon distribution of 
PSR B2224-I-65 in the merged Chandra image is 
slightly deviated from a point source with the de- 
formation also toward the direction behind the 
pulsar proper motion. Although the evidence for 
the extended source cannot be conclusive with the 
existing data, the compact feature suggested by 
the Chandra data can possibly be the portion of 
the highest surface brightness. Complex struc- 
tures of bow-shock PWNe have been observed in 
other systems. For example, XMM-Newton ob- 
servation of a nearby pulsar PSR B1929-I-10 has 
revealed a long X-ray tail of several arc-minutes 
(Becker et al. 2006). Apart from confirming the 
detection of this long trail, a follow-up Chandra 
observation has further resolved a brighter but 
more compact feature with an extent of ~ 10" 
around the pulsar (Hui & Becker 2008). For a 
further investigation of PSR B2224-h65, a dedi- 
cated deep Chandra observation, which simultane- 
ously provides a sub-arcsecond spatial resolution 
and a photon statistic comparable with this XMM- 
Newton observation, is the only way to confirm or 
refute this suggested feature. 

While this compact feature is still questionable, 
the evidence for the ~ 2' long jet is unambigu- 
ous. Although all the X-ray studies suggest its 
non-thermal nature, its physical origin remains to 



be obscure. Very recently, a coUimated ^ 9' long 
X-ray feature has been discovered from a radio- 
quiet 7-ray pulsar PSR J0357-I-3205 (De Luca 
et al. 2011) which is very similar to the jet of 
PSR B2224-t-65. It is interesting to compare the 
observed properties between these two systems. 
The spin-down luminosity of PSR J0357-I-3205, 
E 6 X 10^^ erg/s, is ~ 6 times higher than 
that of PSR B2224-I-65, but it is the lowest among 
all known non-recycled 7— ray pulsars. The low 
column absorption suggests it is a nearby pul- 
sar with a distance of only a few hundred parsec. 
Assuming there is no inclination with respect to 
the sky plane, this suggests the physical extent of 
the trail associated with PSR J0357-F3205 at the 
order of ^ 1 pc which is similar to the case of 
PSR B2224-I-65. In both cases, there is no evi- 
dence of spectral variation along the feature. For 
a distance of 500 pc, the X-ray conversion effi- 
ciency Lx/E of its trail in 0.5 — 10 keV is ~ 10~^. 
Adopting a distance of 1 kpc for PSR B2224-I-65, 
Lx/E of its jet in the same energy band is found 
to be ~ 10~^ which is an order larger than that 
for PSR J0357-I-3205. This might be why the jet 
of PSR B2224+65 can still be detected in spite of 
its larger distance. 

Following De Luca et al. (2011), we esti- 
mate the highest achievable energy of synchrotron- 
radiating particles injected by PSR B2224-I-65. 
First, we assume the highest energy of the elec- 
trons/positrons can be gained is comparable with 
maximum potential drop in the magnetosphere, 

/ . , \l/2 

which is A"I> = (3i?/2cj for an aligned pulsar 

(Goldreich & Juhan 1969). Similar to the case 
of PSR J0357-I-3205, this results in a maximum 
Lorentz factor of 7max ^ 10^. In the presence of 
a magnetic field, these particles will radiate syn- 
chrotron emission with characteristic frequency of 
Vx = ^"^eB Im^c, where B^^q is the local magnetic 
field strength in the emission region in units of 
microgauss. Adopted a typical magnetic field at 
an order of ~ 5 /iG, 7max implies a synchrotron 
photon energy at the order of hv^ ~ 5 keV. 

One difficulty for the synchrotron scenario is 
that it requires the Lorentz factor of the emitting 
particles to achieve the same order of 7max so as to 
produce X-rays at a few keV. On the other hand, 
inverse Compton (IC) scattering between the rel- 
ativistic particles and the relic photons appears to 
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be more viable in producing the observed X-rays. 
We note that the energy density of the rehc pho- 
tons is at the order of ~ 10~^^ erg cm"'^, which 
is comparable with the magnetic energy density 
for a field strength of ~ 5/iG. Therefore the IC 
power cannot be ignored. In the IC scenario, the 
required Lorentz factor to produce ~ 10 keV pho- 
tons is only at the order of 7 ~ 10'* which is much 
easier to achieve. 

The major difficulty for explaining the X-ray jet 
of PSR B2224-H65 in the context of PWN comes 
from its orientation with respect to the pulsar's 
proper motion. The jet protrudes from the pulsar 
and extends toward northwest, which makes an 
angle of ^ 118° from the proper motion direction 
of the pulsar (cf. Fig. [1]) . This peculiar orientation 
makes the bow-shock interpretation questionable 
as this scenario requires the nebular emission lies 
in an opposite direction of the pulsar's velocity 
vector. In view of this difficulty, different expla- 
nations for the nature of the jet association with 
PSR B2224+65 have been proposed, ranging from 
the magnetically confinement of relativistic parti- 
cles leaking from the bow-shock (Bandiera 2008) 
to a picture similar to the AGN jet outflow (John- 
son & Wang 2010). Nevertheless, there is no con- 
sensus on its physical origin yet. 

Since the extended feature of PSR J0357-h3205 
is akin to that of PSR B2224-I-65, a compara- 
tive analysis between these two systems might 
provide us a deeper insight on the nature of 
their jets. While the distance and the proper 
motion direction of PSR B2224-}-65 are well- 
constrained, these are the most important pa- 
rameters to be determined in further investigation 
of PSR J0357+3205. These would enable a more 
reliable estimations of its energetic and most im- 
portantly the orientation of the jet with respect 
to the pulsar's space velocity, which will allow us 
to ascertain if these two systems have a similar 
physical origin. 

5.2. The origin of PSR B2224+65 

PSR B2224+65 is certainly one of the most 
spectacular neutron stars regarding its motion 
through space as its transverse velocity, v = 
865((i/l kpc) km s^^, is among the highest in 
the currently known pulsar population. On the 
other hand, the well measured bow shock with 
its inclination in the plane of the sky suggests a 



negligible radial velocity component (Chatterjee 
& Cordes 2004) which has been confirmed re- 
cently by Tetzlaff et al. (2009) who suggest the 
origin of PSR B2224-F65 to lie within the Cygnus 
0B3 association or a small massive cluster nearby. 
This origin indicates a small radial velocity of 
Vr = — 21I7Q km s^*. The proposed kinematic 
age of Tkin « 0.8 Myr is consistent with a charac- 
teristic pulsar age of Tchar = 1.1 Myr (Hobbs et al. 
2004; Manchester et al. 2005). 

To advance the analysis we presented earlier 
(Tetzlaff et al. 2009), we attempt to identify a 
runaway star which might have been the former 
companion of this pulsar's progenitor star (we re- 
fer to such runaway stars as binary supernova sce- 
nario (BSS) runaways; Blaauw 1961). Therefore, 
we take all runaway star candidates with full kine- 
matics from the runaway star catalogue of Tetzlaff 
et al. (2011) and compare the past flight paths 
of PSR B2224-I-65 and each runaway star vary- 
ing the observables by performing a Monte-Carlo 
simulation to flnd close encounters between the 
pulsar and the runaway star. Simultaneously, we 
calculate the distance to possible parent associa- 
tions/clusters listed in Tetzlaff et al. (2010) (for a 
full description of the procedure we refer to Tet- 
zlaff et al. 2009 and Tetzlaff et al. 2010)0 

For PSR B2224-h65 we take the following data 
for the equatorial coordinates, right ascension a 
and declination S, the distance to the Sun d and 
the proper motion components ^* (/x* — fia cos S) 
and fis (Harrison et al. 1993; Taylor & Cordes 
1993; Yuan et al. 2010); 

a = 336"?469671, S = 65°593217, 
d = 2 ± 1 kpc, 

^* = 144±3masyr"\ fis = 112±3masyr~^ 

(1) 

Here we adopt d = 2 kpc for the calculation as 
published by Taylor & Cordes (1993) and adopt a 
reasonable error of 1 kpc that covers the distance 
inferred by Cordes & Lazio (2002) and Chatterjee 



PSR B2224+65 reaches a distance to the Sun larger than 3 
kpc (the catalogue limit of the current runaway positions) 
if it is older than 1.3 Myr (assuming Vr = km/s). Since 
the characteristic age (^ 1 Myr) provides an upper limit 
on the true age of young (few Myr) pulsars in most cases, 
it is well possible to identify the associated runaway star 
(the former companion, if it exists) among the catalogue 
runaway star candidates. 
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& Cordes (2004). Among 1703 runaway star can- 
didates in the catalogue with fuU kinematic data, 
only nine appear to show possible close encounters 
with PSR B2224+65 (HIP 97359, 98773, 99283, 
99580, 100308, 100409, 101186, 107418, 109393). 
For two of them, HIP 107418 and HIP 109393, sep- 
arations between the runaway star and the pulsar 
are small at rkin ~ 0, i.e. at present; thus we 
exclude these two stars from our analysis. 

Tracing back the seven remaining runaway stars 
and OB associations/ clusters listed in Tetzlaff et 
al. (2010), we find that four stars could have 
been ejected from their parent cluster within the 
last 3 Myr (w 3 • Tchar) (see Table [3]). For them, 
we simultaneously trace back PSR B2224+65, the 
runaway star and the potential parent association 
varying the observables within their confidence in- 
tervals (here we adopt = ± 100 km s^^ for 
PSR B2224+65) to find the closest approach be- 
tween the pulsar and the runaway star (i.e. the 
minimum separation dmm (''kin)) as well as the as- 
sociated time (rkin) and the distance to the respec- 
tive association centre at that time. 

In only one case - HIP 99580 (=V2011 Cyg) 
with Cyg 0B9 being the parent association - the 
pulsar and the runaway star can be located within 
the association boundaries (radius of Cyg 0B9 
w 40 pc) at the same time in the past. The small- 
est separation dmin between PSR B2224-I-65 and 
HIP 99580 that is found after 15 miUion Monte- 
Carlo runs is 0.4 pc. Between 25 pc and 35 pc 
distance from the centre of Cyg 0B9 0, we find a 
peak of small separations dmin between the pulsar 
and the runaway star. In 2540 simulations the dis- 
tance of the pulsar as well as the distance of the 
runaway star to the centre of Cyg 0B9 falls within 
this range. The first part of the distribution of the 
separations dmin between PSR B2224-t-65 and HIP 
99580 of these runs is consistent with the hypoth- 
esis that both stars once were at the same position 
0.65 ± 0.05 Myr in the past (see Figure [S]). 

Although only a small fraction of runs yield 
small separations dmm between the stars as well 
as the centre of Cyg 0B9 the scenario seems likely 
considering the errors on the parameters involved 



We choose a distance range of 10 pc to make it compara- 
ble to the previously suggested neutron star/ runaway pair 
PSR B1929+10/C Ophiuchi, see Hoogerwerf et al. (2001) 
and Tetzlaff et al. (2010) 



(distance d, proper motion components fi* — 

fia cos S and fis, radial velocity Vr of both star4l3). 
The fraction of successful runs is much smaller 
than in the case of the known neutron star / run- 
away pair PSR B1929-I-10/ C Ophiuchi (about one 
order of magnitude; for PSR B1929-K10/ C Oph 
three million Monte- Carlo runs were performed 
approximately twice the number of successful runs 
compared to 15 million here; see Hoogerwerf et 
al. 2001 and Tetzlaff et al. 2010); however, since 
PSR B2224-f65 and HIP 99580 are far more dis- 
tant from the Sun than PSR B1929-hl0 and C Oph 
(the error on the distance is larger by about one to 
two orders of magnitude), this is to be expected. 

In Figure [71 we show the correlations be- 
tween the astrometric parameters which were 
used as input for the 2540 runs for which both, 
PSR B2224-I-65 and the runaway star HIP 99580, 
were situated between 25 pc and 35 pc from the 
centre of Cyg 0B9. For the runaway star only 
little correlation between tt and can be seen 
and the most probable values of all parameters 
fall clearly within the error bars. Whereas for 
the pulsar there is a strong correlation between 
TT and Vr. This is expected since an increase in 
distance (decrease in tt) can be compensated for 
by a larger radial velocity Vr and vice versa. If 
PSR B2224-h65 and the runaway star HIP 99580 
were at the same place 0.65 Myr ago in Cyg 0B9, 
i.e. once were in the same multiple system, the 
present distance of PSR B2224+65 needs to be 
~ 0.8 kpc. Hence, we can confirm previous find- 
ings by Chatterjee & Codes (2004) which con- 
cludes the pulsar have a distance to the Sun of 
^ 1 kpc. 

If this supernova occurred 0.65 Myr in the past, 
it was located at Galactic coordinates I = 77.4 ± 
0.2 deg, b = 2.5tQ[t de£l, and a distance to Earth 

of dEarth = 9681:28 P^- 

If HIP 99580 is a BSS runaway star, it should 
show signs of stellar evolution in the former sys- 
tem such as a large rotational velocity and an 
enhanced helium abundance caused by momen- 
tum and mass transfer during evolution. It was 



The errors on the motion of Cyg OB9 are also taken into 
account but are negligible compared to the stars. 
These coordinates are as viewed at the present time. (i.e. 
Without the correction for the Sun's motion through the 
Galaxy over the last 0.65 Myr.) 
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claimed by Barannikov (1993) that HIP 99580 is a 
spectroscopic binary, however De Becker & Rauw 
(20040 and Turner et al. (2008) could not con- 
firm this. According to Conti & Ebbets (1977), 
the 05 type star HIP 99580 is a rapid rotator with 
vsmi — 270 km s~^. Penny (1996) investigated 
V sin i of 0-type stars and found that only a small 
number of very massive stars show such a large 
tisini. Hence, it is reasonable to assume that it 
gained its large v sin i during the earlier evolution 
in a binary system. Furthermore, BSS runaway 
stars should be blue stragglers, i.e. they appear 
younger, hence bluer, than their parent association 
(see Hoogerwerf et al. 2001 for other examples). 
In Figure|S]we show the colour-magnitude diagram 
of the potential parent association Cyg 0B9 with 
the data adopted from Garmany & Stencel (1992). 
Apparently, the association consists of two major 
groups of different ages, the younger one being ~ 4 
Myr oldEl If HIP 99580 originated from Cyg 0B9, 
it is clearly a blue straggler. 

Assuming an association age of 4 Myr, the pro- 
genitor star of the neutron star must have been 
3.4 Myr old as it exploded. Thus, its mass was 
(67 ± 15) Mq (Tinsley 1980; Maeder & Meynet 
1989; Kodama 1997). For such large masses, 
formation of black holes is expected rather than 
neutron stars (upper limit « 30 Mq, Heger et 
al. 2003). However, due to mass transfer dur- 
ing binary evolution, even more massive stars (> 
50 — 60 Mq) may eventually produce neutron stars 
instead of black holes (Heger et al. 2003; Belczyn- 
ski & Taam 2008). Also, more massive progenitors 
may also produce larger speeds of the neutron star 
(Burrows & Young 2000) which is consistent with 
the extraordinary speed of PSR B2224-h65. 

The authors would like to thank Ralph Neuhauser 



'We note that these authors quote a different radial veloc- 
ity of Vr = -23.1 ± 8.0 km s'^ for HIP 99580 compared 
to —44 ± 4 km/s (Kharchenko et al. 2007) that we used 
here. In this case, the star would still be a runaway star 
as it was identified also from its tangential velocity (Tet- 
zlafF et al. 2011). However, only a tiny fraction of runs 
would then yield small separations between HIP 99580 and 
PSR B2224-I-65 making the BSS scenario unlikely. But 
note that both radial velocities can be correct if the star 
is a spectroscopic binary, then it would be necessary to 
determine the systematic radial velocity. 
Note that Uyaniker et al. (2001) quote an age for 
Cyg OB9 of 8 Myr. 
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Fig. 1.— The 6' x 6' field-of-view of PSR B2224+65 and its jet as observed by XMM-Newton MOSl/2 
cameras {left panel) and PN camera {right panel) in the energy band of 0.5 — 10 keV. The proper motion 
direction and the proper motion corrected position of the pulsar is illustrated by the black cross and white 
arrow respectively. The insets show the 1' x 1' close-up view around PSR B2224-I-65. All images have been 
smooth by a Gaussian with a — S" . 
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Fig. 2. — The X-ray brightness profile in the energy band of 0.5 — 10 keV along the orientation of the jet as 
observed by XMM-Newton MOSl/2 cameras {upper panel) and PN camera (lower panel). The insets show 
the bins used in computing the profiles in the corresponding data sets. Each bin has a size of 25" x 10". 
The average background level and its la deviation are indicated by horizontal lines which were calculated 
by sampling the source-free regions in the 6' x 6' field-of-view of the corresponding data sets as shown in 
Figure [1] 
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Fig. 3. — a. The X-ray image in the energy band 0.5 — 8 keV centered on PSR B2224-I-65 by merging the 
Chandra ACIS data obtained in different epoch. Top is north and left is east. b. X-ray brightness profile 
in 0.5 — 8 keV as sampled from eight bins along the orientation as illustrated in the inset. The dash line 
indicates the expected profile for a point source. 
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Fig. 4. — The energy spectrum of X-ray emission from the position of PSR B2224-I-65 as observed with 
the PN (upper spectrum) and MOSl/2 detectors (lower spectra) and simultaneously fitted to an absorbed 
power-law model {upper panel) and contributions to the fit statistic {lower panel). 
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Fig. 5. — X-ray spectrum of the entire long jet feature associated with PSR B2224-I-65 as observed with 
the PN (upper spectrum) and MOSl/2 detectors (lower spectra) and simultaneously fitted to an absorbed 
power-law model {upper panel) and contributions to the fit statistic {lower panel). 
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Table 1 

x-ray spectral parameters of psr b2224+65 and its jet.*^ 



S^S^t ^ /oT-lO kev' /o.5-ro kev' X' d.O.f 

10^^ cm~^ 10"^"^ ergs cm'^ 10"^^ ergs cm~^ 



Pulsar + Whole jet 



Pulsar 


-^•3-0.7 


9 9+U.2 


2.4+l:« 3.4+i;^ 


42.2 


49 


Whole jet 




1.2 ±0.1 


9-1^2:8 io-ot^;7 










Whole jet - 


- Region 1 ± Region 2 






Whole jet 


9 D+u.y 


1.3 ±0.1 


8.9+^-« 10.0+^» 


86.2 


91 


Region 1 




1.3 ±0.1 


4 4+^-'* 5 0+^-'^ 

^•^-1.5 ^-^-1.6 






Region 2 




1 q+0.2 


■i.ltli 4.6 ±2.0 






Whole jet + Region 1 + Region 2 (with the knot removed) 


Whole jet 


q 9+l.U 
"5-^-0.8 


1.1 ±0.2 


7.7+^-5 


69.2 


79 


Region 1 
Region 2 




1.1 ±0.2 

1 q+0.2 

-'■•>'-o.i 


4.3t2-« 4.7t^? 
4.612-3 4.6t?j 







^Following the method adopted by Johnson & Wang (2010), these parameters are obtained from 
the joint analysis of the spectra of segments of interest with the column densities in the individual 

models tied together. 

''/0.5-10 keV /o.5-\^o keV represent the observed and absorption-corrected fluxes in 0.5—10 keV 
respectively. 
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Table 2 

x-ray spectral properties and the fluxes of the jet-like feature at different epochs. 



Start date 


Obs. ID. 


Tx /o.S-lO keV 

10^1 cm-^ 10" ergs cm-^ 


i^unabs 
JO. 5-10 keV 

IQ-i^ ergs cm-^ s"! 


Chandra Observations 


2000-10-21 
2006-10-06 


755 
7400 


2.7 (fixed) 1.1 ±0.4 8.7+5;'5^ 
2.7 (fixed) 1.0±0.4 8.7t^^o^ 




XMM-Newton Observation 


2009-07-13 


0604420101 


2.7+i;^ 1.2 ±0.1 9.0^^:^ 


10.0+-^ 



20 



Table 3 

Runaway star candidates with possible close encounters to the Guitar pulsar and a 

POSSIBLE kinematic AGE Tkin,* ^ 3 MYR (w 3 ■ Tchar)- THE POSSIBLE PARENT ASSOCIATION AND THE 
associated Tkin,* ARE GIVEN IN COLUMNS 2 AND 3. COLUMNS 4 AND 5 QUOTE THE STAR'S AGE AND 
MASS M* AS INFERRED FROM EVOLUTIONARY MODELS (SEE TeTZLAFF ET AL. 2011). THE LAST 

Column indicates the spectral type as given in Tetzlaff et al. (2011). 



HIP 


parent 


Myr 


r* 
Myr 


Mq 


SpT 


98773 


Cyg 0B7 




11.6 ±0.5 


15.5 ±0.8 


A2V+... 


99580 


Cyg 0B9^ 


1 0+1° 

-■■•"-1.0 


1.5 ±0.5^^ 


44.6 ±5.1 


05e 


100409 


Cyg 0B9 
Cyg 0B7 
Cep 0B6 


o a + lA 

''■D-o.e 

S 4+1-6 

4 0+1-° 


11.1 ±1.1 


15.5 ±2.7 


Bllb 


101186 


Cyg OBg*^ 


1 +0.8 
■^■^-1.2 


3.9 ±0.1 


37.20.1 


09.5Ia 



^Note that Schilbach & Roser (2008) suggested an origin for HIP 
99580 and HIP 101186 in the cluster NGC 6913 2.4 Myr and 2.9 Myr 
ago, respectively. This cluster lies within the Cyg OBI/ Cyg 0B9 
region. 

''This age does not take into account that the star might be a blue 
straggler. As blue straggler, its age can be consistent with 4 Myr 
(cf. Fig. 9). 
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Fig. 6. — a. Distribution of separations dmin of the 2540 runs for which both objects, PSR B2224+65 and 
HIP 99580 were situated between 25 pc and 35 pc from the centre of Cyg 0B9. Drawn as sohd hne is a 
theoretical curve for a distribution of absolute differences between two Gaussian distributed 3D quantities 
with fi — and a = 5.2 pc (see equation 2 in Tetzlaff et al. 2011). b. Distribution of corresponding flight 
times Tkin in the past since the SN. 
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Fig. 7. — Input parameters (parallax tt or distance d, proper motion components /i* — fia cos 6 and fig, 
radial velocity Vr) for the Guitar pulsar and the runaway star HIP 99580 for the 2540 runs for which both 
objects were situated between 25 pc and 35 pc from the centre of Cyg 0B9. The lines mark the currently 
accepted parameter values (solid) and their la error bars (dashed). 
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Fig. 8. — [B — V)q versus Mv diagram of Cyg 0B9 (dots represent member stars according to Garmany 
& Stencel 1992), the potential parent association of the Guitar pulsar and the runaway star HIP 99850 
(star). The solid line represents the 4 Myr isochrone from Marigo et al. (2008) (for solar metallicity; 
http://stev.oapd.inaf.it/cgi-bin/cmd). [B — V)^ and My of HIP 99580 are derived from its spectral type 
according to Schmidt-Kaler (1982). 
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